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We have measured the W boson mass using the D0 detector and a data sample of 82 pb _1 from the 
Fermilab Tevatron collider. This measurement uses W — > ev decays, where the electron is close to a 
boundary of a central electromagnetic calorimeter module. Such 'edge' electrons have not been used 
in any previous D0 analysis, and represent a 14% increase in the W boson sample size. For these 
electrons, new response and resolution parameters are determined, and revised backgrounds and 
underlying event energy flow measurements are made. When the current measurement is combined 
with previous D0 W boson mass measurements, we obtain Mw = 80.483 ± 0.084 GeV. The 8% 
improvement from the previous D0 measurement is primarily due to the improved determination 
of the response parameters for non-edge electrons using the sample of Z bosons with non-edge and 
edge electrons. 



I. INTRODUCTION 

In the past decade, many experimental results have im- 
proved our understanding of the standard model (SM) ||] 
of electroweak interactions as an excellent representation 
of nature at the several hundred GeV scale ||. Dozens 
of measurements have determined the parameters of the 
SM, including, indirectly, the mass of the as-yet unseen 
Higgs boson. The W boson mass measurement plays a 
critical role in constraining the electroweak higher order 
corrections and thus gives a powerful constraint on the 
mechanism for electroweak symmetry breaking. 

Recently, direct high precision measurements of Mw 
have been made by the D0 §, |, || and CDF § collab- 
orations at the Fermilab pp collider, and by the ALEPH 
§, DELPHI @, L3 and OPAL JhJ collaborations at 
the CERN LEP-2 e^e~ collider. The combined result 
of these measurements and preliminary LEP-2 updates 
§ is M w = 80.451 ± 0.033 GeV. The combined indirect 
determination of Mw [|| from measurements of Z boson 
properties at LEP and SLC, taken together with neu- 
trino scattering studies pT| and the measured top quark 
mass U, is M w = 80.373±0.023 GeV, assuming the SM 
P] . The reasonable agreement of direct and indirect mea- 
surements is an indication of the degree of validity of the 
SM. Together with other precision electroweak measure- 
ments, the W boson measurement favors a Higgs boson 
with mass below about 200 GeV. Measurement of Mw 
with improved precision is of great importance, as it will 
enable more stringent tests of the SM, particularly if con- 
fronted with direct measurement of the mass of the Higgs 
boson, or could give an indication of physics beyond the 
standard paradigm. 

The measurements of Mw in the D0 experiment use 
W bosons produced in pp collisions at 1.8 TeV at the 
Fermilab Tevatron collider, with subsequent decay W — > 
ev. The previous measurements are distinguished by 
the location of the electron in a central electromagnetic 
calorimeter (\r) e \ < 1.1) ||] or the end calorimeters 
(1-5 < |?7e| < 2-5) where 77 is the pseudorapidity, 
r\ = — In tan 9/2, and 8 is the polar angle. The measured 
quantity is the ratio Mw/Mz, which is converted to the 
W boson mass using the precision Z boson mass from 
LEP 0. Decays of the Z boson into e + e~ are crucial for 
determining many of the detector response parameters. 



For all previous D0 W boson mass measurements (and 
for other studies of W and Z boson production and de- 
cay) , electrons in the central electromagnetic calorimeter 
were excluded if they were close to the module bound- 
aries in azimuth (<£>). In this paper, we revisit the cen- 
tral electron W boson analysis, adding these hitherto un- 
used electron candidates that appear near the calorime- 

We use a data sample of 82 
1995 run of the Fermilab 



ter module boundaries 13 
pb _1 obtained from the 1994 



collider. 



II. EXPERIMENTAL METHOD AND EVENT 
SELECTION 

A. Detector 

The D0 detector @ for the 1992 - 1995 Fermilab col- 
lider run consists of a tracking region that extends to a 
radius of 75 cm from the beam and contains inner and 
outer drift chambers with a transition radiation detec- 
tor between them. Three uranium/liquid- argon calorime- 
ters outside the tracking detectors are housed in separate 
cryostats: a central calorimeter and two end calorime- 
ters. Each calorimeter has an inner section for detec- 
tion of electromagnetic (EM) particles; these consist of 
twenty-one uranium plates of 3 mm thickness for the cen- 
tral calorimeter or twenty 4 mm thick uranium plates for 
the end calorimeters. The interleaved spaces between 
absorber plates contain signal readout boards and two 
2.3 mm liquid argon gaps. There are four separate EM 
readout sections along the shower development direction. 
The transverse segmentation of the EM calorimeters is 
0.1x0.1 in Arj x A</>, except near the EM shower maxi- 
mum, where the segmentation is 0.05x0.05 in Ary x A</>. 
Subsequent portions of the calorimeter have thicker ura- 
nium or copper/stainless steel absorber plates and are 
used to measure hadronic showers. The first hadronic 
layer is also used to capture any energy escaping the EM 
layers for electrons or photons. The muon detection sys- 
tem outside the calorimeters is not used in this measure- 
ment, except as outlined in Refs. |^, |5| for obtaining 
a muon track sample used to calibrate the drift chamber 
alignment. 

An end view of the central calorimeter is shown in 
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Fig. |l} There are three concentric barrels of modules; 
the innermost consists of thirty- two EM modules, fol- 
lowed by sixteen hadronic modules with 6 mm uranium 
absorber plates, and then sixteen coarse hadronic mod- 
ules with 40 mm copper absorber plates to measure the 
tails of hadronic showers. All previous D0 W boson 
mass analyses using central electrons have imposed cuts 
on the electron impact position in the EM modules that 
define a fiducial region covering the interior 80% in az- 
imuth of each module. Such electrons will be referred 
to in this paper as 'C or 'non-edge' electrons. The re- 
maining central electrons that impact on the two 10% 
azimuthal regions near an EM module edge suffer some 
degradation in identification probability and energy re- 
sponse, but are typically easily recognizable as electrons. 
We will refer to them as 'C or 'edge' electrons. The 
edge region corresponds to about 1.8 cm on cither side 
of the EM module. Those electrons identified in the end 
calorimeters H are labelled 'E'. The end calorimeters 
have a single full azimuth module and consequently have 
no edges. Dielectron pair samples are denoted CC, CC, 
CC, CE, CE, or EE according to the location of the two 
electrons. 

Cryostat Wall 




FIG. 1: End view of the central calorimeter showing the 
arrangement for electromagnetic (EM), fine hadronic (FH) 
and coarse hadronic (CH) modules. The Tevatron Main Ring 
passes through the circular hole near the top of the CH ring. 

The detailed constitution of the EM calorimeter in the 
vicinity of the edges of two modules is shown in Fig. |^. 
The mechanical support structure for the modules is pro- 
vided by thick stainless steel end plates (not shown); the 
end plates of adjacent modules are in contact to form a 
32-fold polygonal arch. The elements of each module are 
contained within a permeable stainless steel skin to al- 
low the flow of liquid argon within the cryostat. Adjacent 
module skins are separated by about 6 mm. The uranium 
absorber plates extend to the skins, so that any electron 
impinging upon the module itself will pass through suf- 



ficient material to make a fully developed EM shower. 
Within the gaps between absorber plates, G10 signal 
boards are etched on both sides to provide the desired 
r] — <p segmentation for readout. The signal boards are 
coated on both sides with resistive epoxy and held at a 
voltage of 2 kV to establish the electric field within which 
ionization drifts to the signal boards. The resistive coat 
is set back from the ends of the board by about 3 mm 
to avoid shorts to the skin. In the region of this setback, 
the electric field fringing causes low ion drift velocity and 
thus reduced signal size, but the shower development is 
essentially normal as the absorber configuration is stan- 
dard. The hadronic calorimeter modules are rotated in 
azimuth so that the edges of EM and hadronic modules 
are not aligned. 




FIG. 2: Construction of central calorimeter EM modules 
in the region near module boundaries. Signal boards have 
the electrode pads for signal collection; readout boards carry 
traces bringing the signals to the module ends. 

The directions of electrons and their impact point on 
the calorimeter are determined [5| using the central 
drift chamber (CDC), located just inside the calorime- 
ter cryostat. This chamber has four azimuthal rings of 
thirty-two modules each. In each module, the drift cell is 
defined with seven axial sense wires and associated field 
shaping wires. The ring 2 and 4 sense wire azimuthal 
locations are offset by one half cell from those of rings 1 
and 3. Half of the sense wires are aligned in azimuth with 
a calorimeter edge and the other half are aligned with 
the center of a calorimeter module. The drift chamber 
z-coordinate parallel to the beam is measured by delay 
lines in close proximity to the inner and outer sense wires 
of each module, using the time difference of arrival at the 
two ends. 
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B. Triggers 

Triggers for the W boson mass analysis, described in 
more detail in Refs. j|, [|, |J, are derived primarily 
from calorimetric information. For the hardware level 
1 trigger, calorimeter signals are ganged into Ar/ x A<p = 
0.2 x 0.2 towers in both EM and hadronic sections. En- 
ergy above a threshold is required for a seed EM tower. 
The hardware refines this to include the maximum trans- 
verse energy tower adjacent to the seed, and requires this 
combination to exceed a fixed threshold. The correspond- 
ing hadronic tower transverse energy must not exceed 
15% of the EM tower energy. The second level trigger 
refines the information in computer processors using a 
more sophisticated clustering algorithm. At level 2, the 
missing transverse energy (fr) components are formed. 
The W boson level 2 trigger requires an EM cluster and 
]pT above a threshold. The Z boson level 2 trigger re- 
quires two EM clusters. In addition, trigger requirements 
are imposed to ensure an inelastic collision, signalled by 
scintillators near the beam lines, and require the event 
to be collected outside times where beam losses are ex- 
pected to occur 0j. For the offline cuts described below, 
the triggers are 100% efficient @, |§. 

C. Data Selection 

The offline data selection cuts are the same as in the 
previous D0 W boson mass analyses. The variables used 
for event selection are: 

• Electron track direction: The track azimuth of a C 
or C electron is determined from the CDC track 
ccntroid and the reconstructed transverse vertex 
position (determined from the drift chamber mea- 
surement of tracks) . We define the axial track cen- 
ter of gravity in the CDC as z tr k- The track pseu- 
dorapidity is then determined from the difference 
between £ t rk and the EM calorimeter cluster center 
of gravity. 

• Distance of the electron from calorimeter module 
edge: The distance along the front face of the EM 
calorimeter module from the module edge is mea- 
sured by the extrapolation of the line from the event 
vertex through the central drift chamber track cen- 
troid. The azimuthal distance from the module 
edge is denoted d e dge- 

• Calorimeter energy location: rjdet is the pseudora- 
pidity of the EM cluster in the calorimeter, mea- 
sured from the center of the detector. The axial 
position of the EM cluster in the EM calorimeter 
is denoted by z c i us . 

• Shower shape: The covariance matrix of energy de- 
posits in forty lateral and longitudinal calorimeter 
subdivisions and the primary vertex z position are 



TABLE I: Offline selection criteria for central and end elec- 
tron candidates. 



Variable 


Central Electron 


End Electron 


l^det | 


< 1.1 


1.5 — 2.5 


Cshape 


< 100 


< 200 


C"trk 


< 5 


< 10 


EMF 


> 0.90 


> 0.90 


,/iso 


< 0.15 


< 0.15 


A 4 




< 4.0 


l^clus | 


< 108 cm 




\z t ik\ 


< 80 cm 





used to define a chisquare-like parameter, £ s hape, 
that measures how closely a given shower resem- 
bles test beam and Monte Carlo EM showers Jig ]. 

• Electron isolation: the calorimeter energies are 
used to define an isolation variable, /; so = (2?f u ii — 
S cor e) /E COTe , where E COTC is the energy in the EM 
calorimeter within 71=0.2 of the electron direction, 
■Efuu is the energy in the full calorimeter within 
K=0A and K = ^A^ + Acj) 2 . 

• Track match significance: a 2 lk = (As /5s) 2 + 
(A£/<5£) 2 measures the quality of the track match, 
where s is the rcj> coordinate and £ is the z coordi- 
nate for the central calorimeter or radial coordinate 
for the end calorimeter. As and A£ are the differ- 
ences between track projection and shower maxi- 
mum coordinates in the EM calorimeter, and 5s 
and 5£ are the corresponding errors jjj [l| . 

• EM fraction: the fraction, EMF, of energy within 
a cluster that is deposited in the EM portion of the 
calorimeter. 

• Electron likelihood: a likelihood variable, A4, based 
upon a combination of EMF, <7t r k, dE/dx in the 
CDC, and&hape @. 

• Kinematic quantities: the transverse momenta of 
electrons, neutrinos, and the W or Z bosons are de- 
noted pT{e),PT{y),PT(W) or Pt(Z). The px(v) is 
determined from the missing transverse energy in 
the event, as discussed below. The invariant mass 
of two electrons is denoted by m ee . 

The requirements for central and end electrons are given 
in Table | 

The selection criteria for the W and Z boson event 
samples are given in Table [n| Non-edge electrons are 
defined as those with d e d ge /d mo d > 0.1, where d mo d is 
the full width of the module in azimuth. Edge electrons 
are required to have d e d ge /d mo d < 0.1. For the Z bo- 
son sample with two electrons in the central calorimeter, 
both are required to have good tracks in the drift cham- 
ber (i.e. passing the <7trk requirement) if either of them is 
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TABLE II: Event selection criteria for W and Z boson sam- 
ples. 



Variable 


W boson sample 


Z boson sample 


pr(e central) 


> 25 GeV 


> 25 GeV 


pr(e end) 




> 30 GeV 


pr(y) 


> 25 GeV 




Pt(W) 


< 15 GeV 




m ee 




60 - 120 GeV 


pvtx | 


< 100 cm 


< 100 cm 



TABLE III: Event sample sizes. 



W boson sample 


No. events 


Z boson sample No. 


events 


C 


27,675 


CC 


2,012 


C 


3,853 


CC 


470 


E 


11,089 


CC 


47 






CE 


1,265 






CE 


154 






EE 


422 



in a central calorimeter edge region; if both are non-edge, 
only one electron is required to have a good track. For Z 
boson samples with one electron in the end calorimeter, 
the end electron must have a good track, while the cen- 
tral electron is required to have a good track only if the 
electron is in the edge region. 

With these selections, we define three W boson sam- 
ples and six Z boson samples, differentiated by whether 
the electrons used are C, C, or E. The numbers of events 
selected in each sample are given in Table III. 



the pt{s) distribution is mainly broadened by the in- 
trinsic pt(W) distribution, and the PtW) distribution is 
smeared by a combination of both effects. The Z boson 
template variable is the invariant mass, m ee . 

The observed quantities used for W boson reconstruc- 
tion are pr(e) and the recoil transverse momentum, 
{It = T.iETiTii, where hi is the unit vector pointing to 
calorimeter cell i, and the sum is over all calorimeter cells 
not included in the electron region. The electron energy 
in the central calorimeter is summed over a Arj x A0 re- 
gion of 0.5 x 0.5 centered on the most energetic calorime- 
ter cell in the cluster. Note that this region spans 2.5 
modules in azimuth, so always contains several module 
edges irrespective of the electron impact point. For the 
end calorimeter, the electron energy sum is performed 
within a cone of radius 20 cm (at shower maximum), 
centered on the electron direction. In both cases energy 
from the EM calorimeter and first section of the hadron 
calorimeter is summed. 

The neutrino transverse momentum in W boson decays 
is taken to be PtW) = ~Pt{c) ~ &t- The components of 
{It in the transverse plane are most conveniently taken 
as M|| = ut • e and u± — ut • (e x z), where e (z) is the 
electron (proton beam) direction. 

The momentum p(ee) = p(e\) + p{&2) and the dielec- 
tron invariant mass define the dielectron system for the 
Z boson sample. The dielectron transverse momentum 
is expressed in components along the inner bisector axis 
fj of the two electrons, and the transverse axis £ perpen- 
dicular to fj. 

The data are compared with each of the templates in 
turn and a likelihood parameter C is calculated. The set 
of likelihood values at differing boson masses and fixed 
width is fitted to find the maximum value, corresponding 
to the best measurement of the mass. Statistical errors 
are determined from the masses at which ln£ decreases 
by one-half unit from this maximum. 



D. Experimental method 



The experimental method used in this work closely re- 
sembles that of previous D0 W boson mass measure- 
ments. We compare distributions from the W and Z 
boson samples with a set of templates of differing mass 
values, prepared using a fast Monte Carlo program that 
simulates vector boson production and decay, and incor- 
porates the smearing of experimentally observed quanti- 
ties using distributions derived from data. The variables 
used for the W boson templates are the transverse mass, 

m T = ^2p T {e)p T {v)[l - cos(0 e - </>„)], (1) 

and the transverse momenta of the electron and neutrino, 
Pt{&) and pt(v). The three distributions depend on a 
common set of detector parameters, but with different 
functional relationships, so that the measurements from 
the three distributions are not fully correlated. As dis- 
cussed in Ref . [|| , the my distribution is affected most by 
the hadronic calorimeter response parameters, whereas 



E. Monte Carlo production and decay model 

The production and decay model is taken to be the 
same as for the earlier measurements || f|, |j| . The Monte 
Carlo production cross section is based upon a perturba- 
tive calculation p7| which depends on the mass, pseu- 
dorapidity, and transverse momentum of the produced 
boson, and is convoluted with the MRST parton distri- 
bution functions @| . We use the mass-dependent Brcit- 
Wigner function pT with measured total width parame- 
ters T\y and T z to represent the line shape of the vector 
bosons. The line shape is modified by the relative parton 
luminosity as a function of boson mass, due to the ef- 
fects of the parton distribution function. The parameter 
/3 in the parton luminosity function C q q = e~P m " /m ee 
is taken from our previous studies ||, |j . 

Vector boson decays are simulated using matrix ele- 
ments which incorporate the appropriate helicity states 
of the quarks in the colliding protons and antiprotons. 
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Radiative decays of the W boson are included in the 
Monte Carlo model (J] based on the calculation of Ref. 
p9| . Decays of the W boson into tv with subsequent 
t — > evv decays are included in the Monte Carlo, prop- 
erly accounting for the r polarization B . 



F. Monte Carlo detector model 

The Monte Carlo detector model employs a set of pa- 
rameters for responses and resolutions taken from the 
data 0]. Here we summarize these parameters and indi- 
cate which are re-evaluated for the edge electron analysis. 

The observed electron energy response is taken to be 
of the form 



E n 



(2) 



The scale factor a that corrects the response relative to 
test beam measurements is determined using fits to the Z 
boson sample; for the C electrons, a — 0.9540 ± 0.0008. 
The energy offset parameter S correcting for effects of 
uninstrumented material before the calorimeter is found 
from fits to the energy asymmetry of the two electrons 
from Z bosons, and from fits to J/ip — > e + e~ and 
7T° — * 77 — > (e + e - )(e + e~) decays. For C electrons, 
5 = — 0.16^o'2i- There is an additional energy correc- 
tion (not shown in Eq. ||) that contains the effects of 
the luminosity-dependent energy depositions within the 
electron window from underlying events, and also cor- 
rects for the effects of noise and zero suppression in the 
readout. This correction is made using observed energy 
depositions in r\ — <f> control regions away from electron 
candidates. We discuss the modification of the energy 
response parametrization for C electrons below. 
The electron energy resolution is taken as 



0_E_ 

E 



E 



n 



(3) 



where © indicates addition in quadrature. The sampling 
term constant s is fixed at the value obtained from test 
beam measurements, and the noise term n is fixed at the 
value obtained from the observed uranium and electron- 
ics noise distributions in the calorimeter. The constant 
term c is fitted from the observed Z boson line shape. 
The parameter values for C electrons [0 are s = 0.135 
(GeV 1 / 2 ), c = 0.0115±° °°^, and n = 0.43 GeV. The res- 
olution parameters are re-evaluated for C electrons be- 
low. 

The transverse energy is obtained from the observed 
energy using Et = EsinO, where the polar angle is ob- 
tained as indicated in Sec. IIC, with the errors taken 
from the measurements of electron tracks in Z boson de- 
cays. 

The efficiency for electron identification depends on the 
amount of recoil energy, uii, along the electron direction. 
We take this efficiency to be constant for uu < uq and 



linearly decreasing with slope so for uu > uq. The param- 
eters of this model for the efficiency are determined by 
superimposing Monte Carlo electrons onto events from 
the W boson signal sample with the electron removed, 
and then subjecting the event to our standard selec- 
tion cuts. For non-edge electrons, uq = 3.85 GeV and 
.so = —0.013 GeV -1 ; these parameters are strongly cor- 
related Q. Since the properties of electrons in the edge 
region are different from those in the non-edge region, we 
re-examine this efficiency below for the C sample. 

The unsmeared recoil transverse energy is taken to be 



u T = -(R Iec q T ) - Att|| p T (e) + a mh m, 



(4) 



where qr is the generated W boson transverse momen- 
tum; i? roc is the response of the calorimeter to recoil 
(mostly hadronic) energy; Aiiii is a luminosity- and in- 
dependent correction for energy flow into the electron 
reconstruction window; a m b is a correction factor that 
adjusts the resolution to fit the data, and is roughly the 
number of additional minimum bias events overlaid on a 
W boson event; and m is the unit vector in the direction 
of the randomly distributed minimum bias event trans- 
verse energy. The response parameter is parametrized as 
Riec = «rcc + A-cc log qr and is measured using the mo- 
mentum balance in the fj (dielectron bisector) direction 
for the Z boson and the recoil system. The Attii parame- 
ter due to recoil energy in the electron window is similar 
to the corresponding correction to the electron energy, 
but is modified to account for readout zero-suppression 
effects. The recoil response is due to energy deposited 
over all the calorimeter, and thus is not expected to be 
modified for the C electron analysis. 

The recoil transverse energy resolution is parameter- 
ized as a Gaussian response with a ICC = s rccy /ur, mod- 
ified by the inclusion of a correction for luminosity- 
dependent event pileup controlled by the a m b parameter 
introduced above. These parameters are fit from the Z 
boson events using the spread of the f] component of the 
momentum balance of the dielectron-recoil system. Since 
the s rec term grows with p v (ee) while the a m b term is in- 
dependent of p v (ee), the two terms can be fit simultane- 
ously. The recoil resolution parameters are not expected 
to differ for the C and C samples. 



III. BACKGROUND DETERMINATION 

As noted above, the W — > tv — > evvv background is 
included in the Monte Carlo simulation. Because of the 
branching ratio suppression and the low electron momen- 
tum, this background is small (1.6% of the W boson sam- 
ple). The remaining estimated backgrounds discussed in 
this section are added to the Monte Carlo event templates 
for comparison with data. 

The second background to the C W boson sample 
arises from Z — > e + e _ events in which one electron is 
misreconstructed or lost. It is taken to be the same as 
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for the C sample, (0.42 ± 0.08)%, since the missing elec- 
tron is as likely to be an edge electron for both C and 
C samples. Small differences in the shape of this back- 
ground in the case where one Z boson electron falls in 
the edge region give negligible modification to the final 
W boson mass determination. 

The third background for the W sample is due to 
QCD multijet events in which a jet is misreconstructed 
as an electron. This background is estimated by selecting 
events with low using a special trigger which is dom- 
inated by QCD jet production. For events with fix < 15 
GeV, we compare the number of events with 'good' and 
'bad' electrons. Good electrons are required to pass all 
standard electron identification cuts, whereas bad elec- 
trons have track match selection cut cr tr k > 5 and require 
Cshape > 100. We assume that the probability for a jet 
to be misidentified as an electron does not depend on 
$Ti and determine it for both C and C samples. The 
ttit distributions for both C and C samples are shown in 
Fig. H Here, and for the pr(e) and pr{v) distributions, 
the C and C samples are statistically indistinguishable; 
the fraction of background events in the non-edge W bo- 
son sample is (1.3 ± 0.2)%, whereas for the edge sample 
it is (1.5 ± 0.2)%. We use the QCD multijet background 
distribution from the C sample [0 for the C analysis. 
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FIG. 3: Comparison of transverse mass distributions for 
background events to W bosons for C (points with error bars) 
and C (solid histogram). The two distributions are normal- 
ized to the same number of events. 

The background for the Z boson sample is composed 
of QCD multijet events with jets misidentified as elec- 
trons. We evaluate this background from the dielectron 
mass distributions with two 'bad' electrons, one in the 
edge region and one in the non-edge region. We find an 



exponentially decreasing shape of the background as a 
function of m ee with a slope parameter of — 0.064 ± 0.022 
GeV -1 for the CC sample, to be compared with a slope 
of -0.038 ± 0.002 GeV- 1 for the CC sample, so we use 
different background shapes for the two samples. The 
fraction of events in the mass region 70 < m ee < 110 
GeV is (3.7 ±3.6)% for the CC sample and (2.2 ± 1.3)% 
for CC. The CE Z boson background is statistically in- 
distinguishable from the CE Z boson sample, so we use 
the background distribution determined in Ref. || for 
the CE Z boson analysis. 

IV. EDGE ELECTRON ENERGY RESPONSE 
AND RESOLUTION 

A. Determination of edge electron response and 
resolution parameters 

The thirty-two central calorimeter modules are about 
18 cm wide in the r<j) direction at the shower maximum. 
Thus the edge regions defined above are about 1.8 cm 
wide. The Moliere radius Am in the composite material 
of the D0 calorimeter is 1.9 cm. Since electrons deposit 
90% of their energy in a circle of radius 1 \m (and about 
70% within 0.5 Am), the choice was made in all previous 
D0 analyses using central electrons to make a fiducial cut 
excluding electrons within the 10% of the module nearest 
the edge. As noted in Section II, we expect that show- 
ers will develop normally over the portion of the central 
calorimeter module edges where energy can be recorded, 
but that the actual energy seen may be degraded. In 
this section we motivate modified edge electron energy 
response and resolution functions, and describe the de- 
termination of the associated parameters. 

A naive modification for the electron energy response 
and resolution parametrization would use the same forms 
(Eqs. H and ||) employed for the non-edge analyses with 
changed values for some of the parameters. Since the pri- 
mary effect expected as the distance, d d gc , of an electron 
from the module edge varies is the loss of some signal, we 
might consider modified values for the parameter a. Fig- 
ure ||(a) shows the result of a fit for the scale factor a 
in a sample of Z boson events in which one electron is 
in a non-edge region, as a function of the position of the 
second electron. A clear reduction in a is observed in 
the edge bin. When the value appropriate for each bin 
in G? c dgc is used in the analysis for the W boson mass, 
we see a significant deviation of My/ in the edge bin, as 
shown in Fig. |](b) . Modifying both a and the parameter 
c in the resolution function does not improve the agree- 
ment for My/ in different regions. We conclude that this 
simple modification of energy response is inadequate. 

Insight into the appropriate modification to the elec- 
tron response and resolution can be gained by compar- 
ing the Z boson mass distributions for the case of both 
electrons in the non-edge region (CC) to that when one 
electron is in the edge region and the other is non-edge 
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(CC). Figure ||(a) shows both distributions (before any 
energy response scaling), normalized to the same peak 
amplitude. The CC distribution agrees well with the CC 
sample at mass values at and above the peak in the mass 
distribution, but exhibits an excess on the low mass side. 
When the CC distribution is subtracted from the CC dis- 
tribution, the result is the broad Gaussian shown in Fig. 
5(b), centered at about 95% of the mass value for the CC 
sample. 

The data suggest a parametrization of edge electron 
response in which there are two components. The first is 
a Gaussian function with the same response and resolu- 
tion parametrizations as for the non-edge electrons, for a 
fraction (1-/) of the events: 



E m 
~E 



E 



n 

E' 



(5) 



(6) 



and the second is a Gaussian with reduced mean and 
larger width to describe the lower energy subset of events. 
Guided by the data, we take the same functional de- 
scription for the response and resolution parameters for 
a fraction / of events: 



E n 



(7) 



The parameters in Eqs. || and ^ denoted without a tilde 
are those from the previous non-edge W boson mass anal- 
ysis j|. Those with the tilde in Eqs. [j|and || are in prin- 
ciple new parameters for the fraction / of edge electrons 
with reduced signal response. 

The modified response is characterized by a reduction 
in the average energy seen for a fraction of the edge elec- 
trons and on average a reduced EMF for edge electrons. 
A potential explanation for the energy reduction as be- 
ing due to electrons that pass through the true crack 
between EM calorimeter modules is not satisfactory. In 
this case the energy missing in the EM section would be 
recovered in the hadronic calorimeter modules giving the 
correct full electron energy. (We note that there is only 
a 14% increase in the number of W boson electrons (c.f. 
Table III) when the azimuthal coverage is increased by 
25% by including the edge region, indicating that some 
electrons in the true intermodular crack are lost from the 
sample.) 

A more plausible hypothesis is that the electrons in 
the edge region shower in the EM calorimeter normally, 
but for the subset of electrons which pass near the mod- 
ule edge, the signal is reduced due to the smaller electric 
drift field in the edge region. In this case too, the average 
EMF is reduced due to the loss of some EM signal, but 
the overall energy is lowered as well. This picture of the 
energy response agrees with the observed behavior seen in 



<Je_ 
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E' 



(8) 
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FIG. 4: Distributions for C samples as a function of the ratio 
of the electron impact distance d e d ge from the module edge 
to the total module width, d mo d: (a) the fitted scale factor a, 
and (b) the fitted W boson mass using the appropriate scale 
factor for each d e d ge bin. The errors are statistical only. 
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FIG. 5: (a) Dielectron mass distributions for CC and 
CC samples, with the CC distribution scaled to give the same 
peak value as for the CC distribution. The solid histogram is 
for the CC Z bosons and the points are for the CC Z bosons, 
(b) The difference between CC and normalized CC samples. 
The curve is a Gaussian fit; no backgrounds are included in 
the fit to the difference. 
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Fig. 5. Our model is probably oversimplified, since even 
within the edge region there can be a range of distances 
between shower centroid and the module edge where the 
electric field is most affected, leading to variable amounts 
of lost signal. The distribution of Fig. 5(b) however in- 
dicates that a single extra Gaussian term in the response 
suffices to explain the data at the present level of statis- 
tical accuracy. We speculate that the convolution over 
impact position contributes to the rather large width of 
the lower energy Gaussian term, relative to that for the 
full energy Gaussian. 

The representation above for edge electron response 
and resolution introduces six potential new parameters: 
a, 5, J, c, n and /. We expect n — n since the electronics 
noise should be unaffected near the edge of a module. 

Since there is no difference in the amount of niate- 
rial before the calorimeter, we would expect that 6 = 6, 
The determination of 6 can be made from the Z boson 
sample data. For the form of the energy response func- 
tion adopted above, the observed Z boson invariant mass, 
m ee , should be 

m ee = ahlz +Fz6 (9) 

in the case that 5 -C E(ei) + E{e%). Here, Mz is the true 
Z boson mass taken from LEP measurements [0 {Mz = 
91.1875 GeV), T z = [E{e l )+E{e 2 )]{l-msw)/m ee , and 
ui is the opening angle between the two electrons e\ and 
€2- Fitting the dependence of m ee on Tz gives 6. We 
find that the Tz dependence for the CC Z boson sample 
is consistent {\ 2 = 8.9 for 9 degrees of freedom) with 
that for the CC Z boson sample, and thus take 6 = 5. 

We argued above that, because the structure of the ab- 
sorber plates extends well past the region where the high 
voltage plane ends, we would expect the same sampling 
constants in edge and non-edge regions. We check this 
hypothesis by dividing the CC Z boson sample into two 
equally populated bins of edge electron energy, E e < 41 
GeV and E e > 41 GeV, for which the mean energies are 
36 and 47 GeV respectively. Using the non-edge value 
of s for both subsamples, we show in Fig. ^ the Z bo- 
son mass distributions and the Monte Carlo expectation 
for the best template fit described in more detail below. 
We find the fitted Z boson masses are 91.10 ± 0.32 GeV 
{E e < 41 GeV) with \ 2 = 4.5 for 14 degrees of freedom 
and 91.06 ± 0.27 GeV {E e > 41 GeV) with x 2 = 12 for 
16 degrees of freedom. The consistency and goodness of 
fit leads us to take s = s. 

We simulate the response of the calorimeter to elec- 
trons in the edge region, using the GEANT [M program 
with all uranium plates and argon gaps included. The 
simulation lacks some details of the actual calorimeter, 
including some of the material between calorimeter mod- 
ules, and contains an incomplete simulation of the de- 
tailed resistive coat pattern on the signal readout boards. 
The resulting distribution of energy for 40 GeV electrons 
impacting upon the edge region of the calorimeter mod- 
ules is shown in Fig. |7]. The Monte Carlo distribution 



closely resembles that seen in the data, with a fraction 
of events showing a broad Gaussian with lower average 
response than the main component of electrons. Within 
the imperfect simulation of calorimeter details, the agree- 
ment with the data is good. The Monte Carlo distribu- 
tion can be well fit with the same functional form (Eqs. 
|8|) used for the data. 

Thus, we conclude that for the C electrons, we must 
introduce only three new parameters a, c and /. In prin- 
ciple, we expect that these parameters may be correlated. 
Our fitting procedure is to first fit the CC Z boson mass 
distribution with uncorrelated free parameters 3, c and 
/. We use the resultant value / = 0.31 as input to a two- 
dimensional binned likelihood fit of the templates to the 
data created by the Monte Carlo, varying both a and c. 
The two-dimensional contours show that the correlation 
between a and c is very small. Thus in the vicinity of 
the maximum likelihood in the two-dimensional fit, we 
can fit one-dimensional distributions for each parameter 
separately. The one-dimensional fits for a and c are re- 
peated iteratively after modifying the other parameter; 
the process converges after one iteration. The results of 
these fits, shown in Fig. |8|, give a = 0.912 ± 0.018 and 
c = O.IOIIq'q^. For these best fit a and c, we make 

a one-dimensional fit for / as shown in Fig. |^ and find 
/ = 0.346 ±0.076. 

To verify that the non-edge scale factor a and the nar- 
row Gaussian width from the non-edge electrons are in- 
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FIG. 6: Z boson mass distributions for (a) edge electrons 
with Et > 41 GeV and (b) edge electrons with Et < 41 
GeV. The histograms are the best fit distributions from the 
Monte Carlo. The curve at the bottom of (b) represents the 
background. 
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FIG. 7: Monte Carlo simulation of the energy response func- 
tion for 40 GeV electrons in the edge region. The points repre- 
sent the Monte Carlo data and a fit using the parametrization 
of Eqs. pi is given by the curve. 



# 6 
4 
2 


= 4 

i 

3 
2 
1 



0.86 0.88 0.9 0.92 0.94 



a 



(51 



0.06 0.08 0.1 0.12 0.14 0.16 0.18 

c 

FIG. 8: (a) Fits to a with edge electron parameters c and 
/ fixed near their optimum values; (b) fits to c with edge 
electron parameters a and / fixed near their optimum values. 
The curves are best-fit parabolas. 



deed appropriate for the fraction (1-/) of edge electrons 
represented with standard response, we perform a fit to 
the CC Z boson sample in which both narrow and wide 
Gaussian parameters are allowed to vary. The resulting 
values for a and oe for the narrow Gaussian are consis- 
tent with those obtained in the non-edge analysis [Q. 

We also look for a dependence of the response param- 
eters on the electron selection variables EMF, /j so , £ s hape 
and Otrk by breaking the Z boson sample into bins of 
each of these variables and fitting for the edge fraction 
/ within each bin. No significant variations are seen. 
The largest is a one-standard-deviation slope in the fit- 
ted / vs EMF distribution, and we examine the effect of 
this small dependence as a cross-check below. 

The resulting likelihood fit to the CC Z boson mass us- 
ing the parametrization given above is shown in Fig. 
For this fit, a set of Z boson events is weighted in turn 
to correspond to templates of Z boson samples spaced at 
10 MeV intervals. The best fit yields M z = 91.20 ± 0.20 
GeV, with a,\ 2 — 10.4 for 19 degrees of freedom. The fit- 
ted Z boson mass agrees very well with the input Z boson 
mass from LEP || used in establishing the parameters 5, 
c and /. The small, statistically insignificant, deviation 
from the input value occurs since we use the values of 
parameters a and c from Ref. [Q and not those which 
give the absolute minimum \ 2 when these parameters are 
varied in the CC analysis. 

We also investigate alternate parametrizations for the 
edge electrons involving a Gaussian-like function with 
energy-dependent width or amplitude. If we adopt the re- 
quirement that such parametrizations add no more than 
three new parameters, as for our choice above, we find 
such alternatives to be inferior in their ability to repre- 
sent the Z boson mass distribution. 



B. Cross checks for edge electron response and 
resolution parameters 

We noted above that the fraction / of reduced response 
electrons in the edge region displays some dependence 
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FIG. 9: Fits to / with edge electron parameters a and c fixed 
at their optimum values. The curve is a best fit parabola. 
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upon the fraction of the total energy seen in the EM sec- 
tion. Thus our fitted parameters have been averaged over 
a range of EMF values. To check that this averaging is 
acceptable, we perform analyses separately on approxi- 
mately equal-sized subsets of events with low and high 
EMF fractions (EMF < 0.99 and EMF > 0.99), for both 
the CC Z and C W boson samples. (Values of EMF 
> 1 are possible due to negative noise fluctuations in the 
hadron calorimeter energy.) For the CC Z boson sample, 
no EMF requirement is made on the C electron. Since 
the values of the Z boson mass in the low and high EMF 
CC Z boson sample subsets differ slightly, and the en- 
ergy scale parameter a for non-edge electrons is used in 
the edge electron response function, we determine the 
appropriate ct's for the two EMF ranges of the CC data 
separately. The relative change for the scale factor a 
for the low EMF non-edge electrons is —0.17%, and for 
the high EMF selection is +0.32%. Using these modi- 
fied values for a, we fit the edge electron parameters 5, 
c and / for each subrange separately. Using these re- 
sults, we create templates using the modified parameters 
and fit for the W and Z boson masses in both subranges. 
The transverse mass distribution was used to obtain My/. 
Table |^ shows the fitted parameters and the resultant 
mass fits for low and high EMF subsets. The W and Z 
masses agree between the two subsets; the difference in 
the fitted Z boson mass between the high and low EMF 
subsets is — 0.47± 0.39 GeV, and for the W boson mass 
is 0.62 ± 0.45 GeV. As expected, the fraction / is larger 



TABLE IV: Fitted parameters for edge electrons, and W and 
Z boson mass values, for separate low and high EMF fraction 
subsamples. 
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FIG. 10: (a) Best fit to the CC Z boson mass distribu- 
tion using the parametrization discussed in the text for edge 
electron response and resolution. The lower curve is the ex- 
pected background, (b) The likelihood function as a function 
of hypothesized Z boson mass. 





Low EMF subset 


High EMF subset 




a 


0.922 ± 0.025 


0.888 ± 0.024 




c 


0.163 ±0.026 


0.047 ± 0.027 
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0.45 ± 0.08 


0.25 ±0.06 


M w 


(GeV) 


80.23 ± 0.34 


80.84 ± 0.29 


M z 


(GeV) 


91.43 ± 0.31 


90.96 ± 0.25 



for the low EMF subset, and the width parameter of the 
Gaussian resolution c is larger. The errors quoted are sta- 
tistical only; we estimate that inclusion of the systematic 
errors would roughly double the total error. We conclude 
that the analyses for the two subsets in EMF are in good 
agreement, validating our choice to sum the two samples 
in the primary analysis. 

The averaging over the range of EMF values that oc- 
curs in our analysis is acceptable if the electron EMF 
distribution is the same for the C W boson sample and 
the Z boson CC sample used to obtain the parameter 
values. Fig. |ll| shows the EMF distributions for these 
two samples overlaid; they are statistically consistent. 
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FIG. 11: EM fraction distribution of edge electrons for the 
C W boson (data points) and CC Z boson (histogram) sam- 
ples. The Z boson sample is normalized to the W boson 
sample. 

The parameters for edge electrons discussed above are 
determined from the CC Z boson sample. It is thus use- 
ful to examine other samples in which C electrons par- 
ticipate to demonstrate the validity of the parametriza- 
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tion. The CE dielectron sample with one edge central 
calorimeter electron and one end calorimeter electron, 
using the energy response and resolution of Ref. || for 
the end electrons, is shown in Fig. OL This distribution 
is fit with Z boson mass templates and yields the result 
M z = 91.10±0.42 GeV (statistical) with X 2 = 9.8 for 13 
degrees of freedom, in good agreement with the precision 
LEP Z boson mass determination. When the reduced 
response term for a fraction / of central electrons in the 
edge region is omitted, the fitted Z boson mass is about 
one standard deviation low, and the quality of the fit 
deteriorates to x 2 — H-7. 

We also examine the dielectron sample in which both 
electrons are in the central calorimeter edge region. The 
data shown in Fig. |l3| comprising 47 events is fitted to 
Z boson mass templates to give Mz = 90.38 ± 0.33 GeV 
(statistical). The fit gives \ 2 = 8-5 for 6 degrees of free- 
dom. When the systematic errors are included, this result 
is in reasonable agreement with the LEP precision value 
for M z - 
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FIG. 12: Best fit to the CE Z boson mass distribution us- 
ing the parametrization discussed in the text for the central 
calorimeter edge electron response and resolution, and the 
parameterization of Ref. || for the end calorimeter electron. 
The histogram is the best fit from Monte Carlo, and the lower 
curve is the background. 

As a final cross check, we subdivide the full Z bo- 
son sample into five subsets, in which one electron (the 
'tagged' electron) is required to be in a bin determined 
by the distance d e dge from the nearest module edge. Five 
equal-sized bins span the range < ii e dge/dmod < 0.5. 
The other electron is required to be in any of the non- 
edge bins not populated by the tagged electron. A com- 
panion sample of W boson candidates, subdivided into 
the five dedge bins, is also formed. For each of the Z 



boson samples, the tagged electron response is fitted as 
described above with a variable energy scale factor a us- 
ing the LEP precision value as input. This modified scale 
factor is then used for the W boson subsamples to ob- 
tain a best fit W boson mass. The results are shown in 
Fig. |l4], where the points in the bin < dodgo/^mod < 0.1 
are those from the edge electron with additional param- 
eters as described above. The resulting W boson mass 
values are consistent over the five bins, indicating that 
our energy response correction analysis is acceptable. 



V. OTHER PARAMETER DETERMINATIONS 

Although we expect that the main modifications to the 
previous non-edge electron W boson analyses are the re- 
sponse and resolution parametrizations discussed in Sec- 
tion [Fv|, there are some other parameters that could be 
sensitive to the location of the electron relative to the 
module boundary. 

The observed electron and recoil system energies are 
changed from the true values by the energy from the un- 
derlying event deposited in the region used to define the 
electron. This component of energy must be subtracted 
from the observed electron energy and added to the re- 
coil. In Ref. |4| we found this correction to be dependent 
on the electron rapidity and on the instantaneous lumi- 
nosity. The size of the region used to collect the electron 
energy is A77 x Acj) = 0.5 x 0.5, spanning two and a half 
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FIG. 13: Best fit to the CC Z boson mass distribution us- 
ing the parametrization discussed in the text for the central 
calorimeter edge electron response and resolution. The his- 
togram is the best fit from Monte Carlo, and the lower curve 
is the background. 
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times the size of a module in the <p direction. Thus the 
underlying event correction can only be very weakly de- 
pendent on the location of the center of this region, and 
we take the correction to be the same as for the non-edge 
analysis. Also, the recoil system has its momentum vec- 
tor pointing anywhere in the detector in both the edge 
and non-edge analyses. Thus we do not modify the pre- 
vious parameters controlling the recoil system response 
and resolution. 

The efficiency for finding electrons is modified by the 
underlying event energy within the electron region. The 
efficiency depends on u\\ , since when there is substantial 
recoil energy near the electron, the isolation requirement 
will exclude more events than when the recoil energy is 
directed away from the electron. Since the electron en- 
ergy itself is modified near the module edge, this effi- 
ciency could be different for C and C electrons. To inves- 
tigate this effect, we compute the average /i so for both 
C and C samples. We find that (/i so ) for the C sample 
is 1.08 ± 0.15 times that for the C sample. We expect 
about a 3% increase in (/i so ) since its definition involves 
the EM energy near the core of the shower, which is re- 
duced for C electrons. A modified distribution of /j so 
can only affect the U\\ efficiency if there is a change in 

the u\\ distribution in the C events relative to that for 
the C electrons. We see no difference in the (/i so ) value 
in hemispheres Km < and un > for the C events. This 
observation, and the statistically insignificant difference 
for (/i S o) for C and C samples, lead us to retain the pre- 
vious parametrization for the u» efficiency. 
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FIG. 14: (a) The EM scale factor a and (b) the fitted W 
boson mass in bins of d c d go /d m od using response parameters 
from a Z boson sample requiring one electron in the same bin. 



Since photons radiated from electrons are found domi- 
nantly near the electron, these photons also populate re- 
duced response regions in the edge electron analysis. For 
our analysis we have chosen to generate such radiation 
with the response parameters found for the C electrons. 
However some of the radiated 7's strike the non-edge 
region and should thus be corrected with the non-edge 
response. We calculate that the difference between the 
photon energy using the edge response and a properly 
weighted response across the module is only 3.5 MeV, 
resulting in a negligible less shift in the W boson mass 

ill- 

When an electron impacts the calorimeter near a cell 
boundary, as occurs near the module edge, its position 
resolution in r<f> is improved typically by about 20% [|l4| . 
This means that the determination of the electron clus- 
ter azimuth is more accurate for C than for C electrons. 
The effect of improved azimuthal precision in the C sam- 
ple has however been incorporated by fitting the energy 
response and resolution parameters for the CC Z boson 
sample, so no additional correction is needed. 

The small modification to the electron energy (a 4% re- 
duction in 35% of the electrons in the edge region) could 
affect the trigger efficiency near the threshold. We deter- 
mine that this effect is negligible. 



VI. W BOSON MASS DETERMINATION 

A. Mass fits 

Monte Carlo templates are prepared for the W bo- 
son transverse mass tot, electron transverse momentum 
Pr{e), and neutrino transverse momentum pt{v), using 
the production, decay, and detector parameters discussed 
in Sections II and IV. The estimated backgrounds de- 
scribed in Section III are added to the Monte Carlo W 
boson decays. Families of templates are made for W 



boson masses varied in 10 MeV steps between 79.6 and 
81.6 GeV. The templates are compared to the data in 
the ranges 60 < m T < 90 GeV, 30 < p T (e) < 50 GeV, 
and 30 < p T {v) < 50 GeV, with bins of 100 MeV for 
transverse mass and 50 MeV for the transverse momen- 
tum distributions. For each specific template with fixed 
Mw, we normalize the distributions to the data within 
the fit interval and compute a binned likelihood 



N 



C(m) = JJp"*(to) 



(10) 



where Pi(m) is the probability density for bin i with the 
W boson mass taken as to, rii is the number of data 
events in bin i, and N is the number of bins in the fit 
interval. We fit — ln£(m) with a quadratic function of 
to. The value of m at which the function assumes its 
minimum is the fitted value of the W boson mass and 
the 68% confidence level statistical error corresponds to 
the interval in m for which — In/! (to) is within half a unit 
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TABLE V: Fitted W boson masses and x 2 /degrees of free- 
dom. 



Distribution 


Fitted mass 


X 2 /d.o.f. 


rriT 


80.596 ± 0.234 


45/29 


Pr{e) 


80.733 ± 0.263 


38/39 


Pt{v) 


80.511 ±0.311 


45/39 



of the minimum. The best fit tot, Pt{&) and pr(y) dis- 
tributions and the associated likelihood curves are shown 
in Figs, [l^ - [l?]. The fitted values for Myy and x 2 from 
each of the distributions are given in Table 0. The errors 
shown are statistical only; the values of My/ obtained 
from the three distributions are in good agreement. 

We study the sensitivity of the fits to the choice of 
fitting window by varying the upper and lower window 
edges by ±10 GeV for the transverse mass and by ±5 
GeV for the transverse momentum fits. Figure [l8] shows 
the change in My/ as the upper and lower window edges 
for the transverse mass fit are varied. The shaded bands 
correspond to the 68% probability contours, determined 
from an ensemble of Monte Carlo W boson samples with 
the chosen window edges. The dashed lines indicate the 
statistical error for the nominal fit. The points for dif- 
ferent window edges are correlated, as the data with a 
larger window contains all the data in a smaller window. 
The deviations of My/ are in good agreement for differ- 
ing choices of window. Similar good agreement is seen in 
varying the windows for the pr(e) and pr(y) fits. 



B. Mass error determination 

In addition to the statistical errors determined from 
the fits, there are systematic errors arising from the un- 
certainties in all of the parameters that enter in the 
Monte Carlo production, decay and detector model. 
These parameters, summarized in Table VI, form a pa- 
rameter vector P. The definition and determination of 
the parameters are described above and in Ref . Q . The 
recoil response takes into account the joint effects of two 
correlated parameters a rcc and /3 rcc - We assign an uncer- 
tainty in My/ for the uncorrelated errors obtained from 
the principal axes of the a roc — /3 TO c error ellipse Q . The 
recoil resolution depends on correlated parameters s rcc 
and a m b (§, an d the U|| efficiency depends on corre- 
lated parameters uq and so; these correlated pairs are 
treated similarly to those for the recoil response. The 
set of production model errors include the parameters 
due to parton distribution function (PDF) uncertainty, 
W boson width plj , the parameters determining the W 
boson production px spectrum, and the parton luminos- 
ity function. We take the components of the production 
model error to be uncorrelated. The PDF error is taken 
from the deviation of the W boson mass comparing || 
MRS (A)' p§, MRSR2 f2§, CTEQ5M §1, CTEQ4M 



f2j§, and CTEQ3M [g§ PDF's to our standard choice of 
MRST. 

In all, we have identified Np = 21 parameters that 
determine the model for the Monte Carlo: the eighteen 
used in the previous studies and the three new parame- 
ters related to the edge electrons (5, c and /). 

The parameters Pi are determined from Ny = 32 aux- 
iliary measurements using several data sets which in- 
clude the CC and CC Z boson samples, special mini- 
mum bias and muon samples for determining drift cham- 
ber scales and underlying event properties, and external 
data sets that are used to constrain the W boson produc- 
tion model. The measurements using these special data 
sets are denoted Y] (I — 1, ... Ny) with uncertainties aj. 
Each measurement puts constraints on one or more of the 
parameters Pi. Measurements Yj are related to the pa- 
rameters Pi through the functional relation Yj = Fj(P). 

We form the x 2 for the set of measurements 



E 

i,j=i 



[Yj-FjiP^CYjr^Yj-FjjP)} 



KM) 



(ii) 



where Cjj = (AYjAYj) is the covariance matrix of the 
measurements, determined from Monte Carlo calcula- 
tions. If the deviations of the measurements from their 
means are taken to be linearly related to the parameters 
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FIG. 15: (a) Comparison of the data (points) and the Monte 
Carlo predicted distribution (histogram) in transverse mass, 
using the fitted value for Mw- The Monte Carlo distribution 
is normalized in area to the number of W boson events within 
the fitting window. The estimated backgrounds are indicated 
by the lower curve, (b) The distribution of calculated likeli- 
hood values as a function of the assumed W boson mass. The 
curve is a fitted parabola. 
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in the region of the \ 2 minimum: 



AY, 



Np 

E 

3=1 



(12) 



where Djj = dF I /dP 3 , the minimum of the \ 2 can be 
found analytically. The parameter covariance matrix Cf^ 
can be then calculated from Cjj and the derivatives Djj . 

This analysis is carried out for the three distinct mea- 
surements of Mw for the edge electrons (tot, PT(e) and 
PTiy)). Each measurement depends on the set of param- 
eters, P, discussed above. For the Nm = 3 separate mass 
measurements m a (a = 1, ...Nm), the mass measurement 
covariance matrix C^p is obtained from 



Np 

E 

k,l=l 



n M r P n M 



(13) 



where = dm a /dPj. The correlation of the statistical 
errors is obtained from studies of Monte Carlo ensembles; 
these correlations are shown in Table 



VII 



We can fit for the best combined mass value Mw by 
minimizing the \ 2 27 



N M 

x 2 - E 

a, 0=1 



(m Q - M w )H af3 (m - M w ), (14) 
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FIG. 16: (a) Comparison of the data (points) and the Monte 
Carlo predicted distribution (histogram) in electron trans- 
verse momentum, using the fitted value for Mw- The Monte 
Carlo distribution is normalized in area to the number of W 
boson events within the fitting window. The estimated back- 
grounds are indicated by the lower curve, (b) The distribution 
of calculated likelihood values as a function of the assumed 
W boson mass. The curve is a fitted parabola. 
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FIG. 17: (a) Comparison of the data (points) and the Monte 
Carlo predicted distribution (histogram) in neutrino trans- 
verse momentum, using the fitted value for Mw The Monte 
Carlo distribution is normalized in area to the number of W 
boson events within the fitting window. The estimated back- 
grounds are indicated by the lower curve, (b) The distribution 
of calculated likelihood values as a function of the assumed 
W boson mass. The curve is a fitted parabola. 
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FIG. 18: Variation of the fitted W boson mass with (a) 
the lower edge and (b) the upper edge of the fit window for 
the transverse mass distribution. The shaded regions and the 
dashed lines are described in the text. 
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TA "RT "FT VT 


Parameters P used in the W boson mass deter- 


mination 




Parameter 


Description 


a 


EM energy response scale for non-edge e 


a 


EM energy response scale for edge e 


S 


EM response offset 


c 


EM resolution constant for non-edge e 


c 


EM resolution constant for edge e 


7 


fraction of low response e in edge region 


Pcdc 


drift chamber position scale factor 




recoil energy response scale constant 


f3rec 


recoil energy response scale Q 2 dependence 


Srec 


recoil energy resolution 


Otmb 


recoil energy from added minimum bias events 


AU|| 


underlying event energy correction in e window 


uo 


U|| cutoff for constant efficiency 


so 


slope of u\\ efficiency vs u\\ 


&w 


background to W boson distribution 


r 7 


coalescing radius for photon radiation 


2^ 


fi'TTCiv for 2 *v radiation 


PDF 


error from varying PDF 


r»- 


W boson width 


P 


parton luminosity 


92 


Q 2 dependence of W boson production 



TABLE VII: The statistical correlation coefficients for the 
three measurements of the W boson mass. 

rriT Pr{e) Pt(v) 

m T 1 0.669 0.630 

p T (e) 0.669 1 0.180 

p T {v) 0.630 0.180 1 



TABLE VIII: Errors (in MeV) for the three W boson mea- 
surements 



Source 


TTbT 


TjTi [PI 




tit Q"f"l 


234 


263 


311 


ppIctp T^IVT qppIp ( rv i 


265 


309 


346 


CO EM scale (a) 


128 


131 


113 


CC EM offset (A\ 


142 


139 


145 


palnrirnpfpr iinifnrmitv 


10 


10 


10 


CVlC scale 


38 


40 


52 


r~\Ct c\r nTni inno 
UeLLxVtLI U LillUo 


1 n 


on 


on 


CC EM constant term c 


15 


18 


2 


edge EM constant term (c) 


268 


344 


404 


TT*tl ct 1 OjTI PiT P^7"PTli"<5 1 f 1 
HdLllUll Ul CVClluo I J 1 


g 


14 


22 


hadronic response 


20 


16 


46 


hadronic resolution 


25 


10 


90 


uii correction 


15 


15 


20 


u\\ efficiency 


2 


9 


20 


parton luminosity 


9 


11 


9 


radiative corrections 


3 


6 


< 1 


27 


3 


6 


< 1 


Pt{W) spectrum 


10 


50 


25 


W boson width 


10 


10 


10 



TABLE IX: The full correlation coefficients for the three 
measurements of the W boson mass. 





mT 


p T {e) 


Priy) 


rriT 


1 


0.90 


0.89 


p-r(e) 


0.90 


1 


0.76 


p T {v) 


0.89 


0.76 


1 



where H = (C ) . The best fit is given by 

N M N m 



a, (3=1 



a, /3=1 



with error 



TV, 



a, 13=1 



(16) 



The resultant W boson mass measurements using elec- 
trons in the edge region are 

M w = 80.596 ± 0.234 ± 0.370 GeV (17) 
for the mr(W) fit, 

M w — 80.733 ± 0.263 ± 0.460 GeV (18) 
for the Pt{&) fit, and 

M w = 80.511 ± 0.311 ± 0.523 GeV (19) 



for the pt{v) fit, where the first error is statistical and 
the second is systematic. The breakdown of the contri- 



butions to the systematic errors is shown in Table VIII 



The PDF error is taken as the difference on the com- 
bined W boson mass between the CTEQ3M and MRST 
choices, for which my/ differs maximally. The com bined 
mass error from this source (not shown in Table VIII ) 
is 19 MeV. The errors associated with the broad Gaus- 
sian parameters in the edge electron response (5 and c) 
dominate the systematic errors. 

The three measurements of Mw are correlated as 



shown in Table IX; when combined taking these correla- 



tions into account, we obtain 



M w = 80.574 ± 0.405 GeV, 



with x 2 = 0.61 for two degrees of freedom. 



(20) 
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VII. COMBINATION OF ALL D0 W BOSON 
MASS MEASUREMENTS 

The analysis presented here for the edge electrons 
brings two new ingredients to the D0 W mass measure- 
ments. First, the edge electron sample is statistically 
independent of all other measurements, and thus can be 
combined to give an improved Mw measurement. Sec- 
ond, the added statistics of the CC and CE Z boson sam- 
ples can be used to refine the knowledge of the electron 
response parameters for non-edge central calorimeter or 
end calorimeter electrons. The improved energy scale 
factors in turn give improved W boson mass precision. 



A. Modified non-edge electron W boson mass 



to be compared with the published value of Mw = 
80.446±0.108 GeV f§. The new end calorimeter electron 
scale factor gives a modified W boson mass: 



M w = 80.679 ± 0.209 GeV, 



(24) 



to be compared with the published value from the end 
calorimeters of M w = 80.691 ± 0.227 GeV @. 

With the modified scale factors for C and E electrons, 
we obtain 



M w = 80.481 ± 0.085 GeV, 



(25) 



with x 2 = 5.5 (6 degrees of freedom) for all non-edge cen- 
tral and end calorimeter measurements, compared with 
the previous determination M w = 80.482 ± 0.091 GeV 



Using the CC sample and the same fitting procedure 
described in Section [V for the C electrons, we have 
obtained a scale factor a = 0.9552 ± 0.0023 for the 
non-edge electrons. This value can be compared with 
the previous determination from the CC sample j| of 
a = 0.9540 ± 0.0008. The correlation matrix for CC and 
CC measurements is calculated in the manner discussed 
in Section VI. 

Similarly, the CE sample can be used to constrain the 
scale factor a for both end and non-edge central elec- 
trons (recall that the central edge electrons contain a 
fraction (1 — / G dgc) of events whose scale factor and reso- 
lution are identical to those of the central non-edge elec- 
trons). Taking into account the correlations, we obtain 
a = 0.9559 ± 0.0107 for electrons in the non-edge region 
of the central calorimeter and a ~ 0.9539 ± 0.0085 for 
the electrons in the end calorimeter. The latter value 
can be compared with the previous value || of the end 
calorimeter electron scale of 0.9518 ± 0.0019. 

Taking the two new measurements of a for the cen- 
tral calorimeter together with the previously determined 
value, we obtain 



a = 0.9541 ± 0.00075. 



(21) 



This new scale factor is higher than the previous value 
by 0.0001, and the error is reduced by 6%. For the end 
calorimeter, the new combined scale factor is 



a = 0.9519 ±0.0018, 



(22) 



again higher than the previous value by 0.0001 with a 5% 
reduction in error. 

In principle, the added data could also improve the 
precision for the resolution constant term c in the central 
and end calorimeters, but in practice it does not. 

With the new values for the scale factors for the non- 
edge central calorimeter electrons, we obtain modified re- 
sults for the non-edge central calorimeter W boson mass: 



B. Combined W boson mass from all D0 
measurements 

With the edge electron mass determinations reported 
in this paper, there are now ten separate D0 W boson 
measurements: the Run la central calorimeter transverse 
mass measurement [^j, three Run lb central calorime- 
ter non-edge measurements Q| (from the transverse mass 
and electron and neutrino transverse momenta), three 
Run lb end calorimeter measurements ||, and the three 
present measurements of the central calorimeter edge 
electrons. Combining these ten mass measurements us- 
ing the method outlined in Sec. VI and an expanded set 



of measurements and parameters to incorporate also the 
end calorimeter electrons, we obtain a final D0 combined 
measured value for the W boson mass of 



M w = 80.483 ± 0.084 GeV 



(26) 



with x 2 — 6-3 (9 degrees of freedom). This value is to be 
compared with our previous || combined measurement 
of M w = 80.482 ± 0.091 GeV. The edge electrons in the 
central calorimeter have improved the precision over the 
previously published results by 7 MeV, or 8%. 



VIII. SUMMARY 

Using a sample of electrons which impact upon the 10% 
of a central calorimeter module closest to either module 
edge in azimuth, we have made a new measurement of 
the W boson mass, and have refined our knowledge of 
the energy scale for previously used electrons that are in 
the interior 80% of the central calorimeter modules or are 
in the end calorimeters. Adding the new measurement 
using the edge electrons gives the final combined result 



M 



w 



80.483 ±0.084 GeV (D0). 



Ah 



w 



80.438 ±0.107 GeV, 



(23) 



Combining the new D0 W boson mass value reported 
here with the CDF || and UA2 p3] measurements, tak- 
ing into account the updated correlated systematic errors 
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for the three experiments due to parton distribution func- 
tion uncertainties and multiple photon radiation gives 




M w = 80.454 ± 0.059 GeV (pp). 

This is an improvement over the previous measurement 
from hadron colliders of M w = 80.452 ± 0.062 GeV @. 
Further combining with the LEP experiments' prelimi- 
nary measurement M w = 80.450 ± 0.039 GeV J|, we 
find the world average W boson mass from direct mea- 
surements to be pjj 

M w = 80.451 ± 0.033 GeV (world). 

The edge electrons used in this analysis represent a 
14% increase in the central calorimeter W boson sample, 
and an 18% increase in the total Z boson sample. The 
larger sample sizes should be of use for all subsequent 
studies of vector bosons in D0. 
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